Analysis of Packing Function 
Solutions for Monomeric 
Proteins 
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START: 



Set Up 
Parallelization 



Initialize 
. Constants and 
Sine/Cosine Tables 



Input 1 & 
Initialize the List of 
Fourier Indices (hid) 



Get Biggest Sphere 

Without Overlap 
of Symmetry Mates 



Input 2 



Confirm 
Calculation Mode 



Initialize 
Tables of Spherical 
Bessel Functions and 
Bessel Function O's 



Determine Limits 
on 1 & m indices and 
on n indices ol each I 



For modes 5 thru 7: 
Read expected SHSB 
coefficient 8c 
deviation values 



For Use with a Radial Correction or 
with Modes 5 thru 7: 



Get 

Fractionalizaiion 
Matrix 



Partition, the List 
of Fourier Indices 
Between 
Proccssor(s) 



Modes 1& 2 (cont'd) 




Cycle 




F 


r Each Choice (mstop) of 


Stopping Value for the SHSB 




m index 






Update Correlation 






Exponentiation 






Factor 





Convert Each Fourier 
Index (hkl) to Polar 
Coordinates & Get 
Its Bessel Argument 



Get the FT of a 
Crystalline Unit Cell 
Filled with Symmetry 
Related Spheres 



Calculate a Radial 
Correction Factor 
for Each *'n M Index 



Calculation Mode-Specific Routines: 

Modes 1 & 2 
(Unplutsed Diffraction Amplitudes to 
Phased FT of SHSB-modeled Unit Cell) 



Mode-Specific 
Input 



If there Is a Prior Model 
(On cycles after the 1st cycle 
of the 1st run of the program): 



On Cycles After the 
1st value of mstop 

Update the Hie # for 
the Calculated FT 



For Each input SHSB Model 



Convert the File # 
& the Model # 
to a File Name 



Read Input 
FT & SHSB coefs 
of the Prior Model 



If there Is No Prior Mode?; 



Model the Crystal as 
a Crystalline Unit Cell 
Filled with Symmetry 
| Related Sphe res j 



Get Input File * for 
Calculated FT, if 
there is a Prior Model 



Get a Starting Value 
(msus) for the SHSB 
in Index for the Next 
Packet of m Values 



Initialize Correlation 
Exponentiation 
Factor 



Update the Model 
by Adjusting the 
Contributions of 
each SHSB function 
for each lmn index in 
the m-packet from 
'msus' to 'mstop* 



FtowQ^lbrlteMBinDihfcil^gnBnfcr-fii^r-: Optima to compute a the FT of a SHSB Model of Crystal 



l.f? 





li 

Figure ~ A schematic example: Two choices for filling the same 
portion of a crystal unit cell from an orthorhombic Spacegroup. 
Although the spheres on the right are smaller than those on the left, 
for some crystals, the local maximum in the packing on the right 
wold be the packing of maximal consistency with the 
crystallographic data. 



Initialize 
F nationalization 
Matrix 



Initialize the Equal 
Partitioning of the 
Fourier (hkl) Index 
between Processors 



On 1st Cycle of 1st Run : 



Prescale Observed 
Diffraction to that 
of a Unit Cell of 
Spheres 



Define the First 
SHSB Index Tripler 
(Imn) for which to 
Consider Model IFI's 



Initialize for Index- 
by-lndex Update 
of Origin-Centered 
SHSB Basis Function 



Modes 4 & 5 only: 
Initialize Buffers for 
Cumulative Update of 
Fourier Representation 



Initialize Pointers to 

Stored Fourier 
Representations of 
Model and of Basis 



Mode 3 only: 
Get File Name from 
File # & Open It to Let 
SHSB Coefs.be Read 



For each "m" Index 
(0 to maximum "m" ) 



For Each hkl in this 

Partition: 
Update "m" Recursion 

Formula for Fourier 
Representation of the 
Origin-Centered SHSB 



For each "I" Index 
(present "m* to maximum T) 



For Each hkl in this 
Partition: 
Update "l,m" Recursion 

Formula for Fourier 
Representation of the 
Origin-Centered SHSB 



For each "n* Index 
(1 to maximum "n" for each *!") 



For Each hkl in this 

Partition: 
Update "n" Recursion 

Formula for Fourier 
Representation of the 
Origin-Centered SHSB 



Depending on Mode: 
Choose the # of Passes 

and # of Presumed 
Phase Angles Needed 
for the SHSB coef. with 
this SHSB index (Imn) 



Set the Presumed 
Amplitude of the 
Origin-Centered 
SHSB Basis Function 



FIRST PASS: 



Initialize Registers: 
Overall Comparison 
of Correlation Coef. 
& Other Statistics 



Renrtialize Pointers to 
Storage Sites for Fourier 

Representations of the 
Full-Unit-Cell SHSB Basis 



Parallel Processor Version: 



Set # of calculations to: 
(# of presumed values of 
SHSB coef.'s phase) 
X 

{P of stored accumulated 
SHSB models for trial 
combination with this 
new SHSB component) 



Given: # of processors 

# of hkl partitions 

# of calculations 
Get: # of required 

rounds of trial 
combinations 



For each round of trial 
combination on this processor 



Single Processor Version: 
(Outer Loop) 



For each presumed value 
of the SHSB coef.'s phase 



Initialize Registers: 
Angular Comparison 
of Correlation Coef. 
& Other Statistics 



